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td Abstract. Realising controlled quantum dynamics via the magnetic interactions 

^ between colour centers in diamond remains a challenge despite recent demonstrations 

, ^, ' for nanometer separated pairs. Here we propose to use the intrinsic acoustical phonons 

in diamond as a data bus for accomplishing this task. We show that for nanodiamonds 
the electron-phonon coupling can take significant values that together with mode 
frequencies in the THz range, can serve as a resource for conditional gate operations. 
Q^ Based on these results we analyze how to use this phonon-induced interaction for 

constructing quantum gates among the electron-spin triplet ground states, introducing 
the phonon dependence via Raman transitions. Combined with decoupling pulses this 
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"^" offers the possibility for creating entangled states within nanodiamonds on the scale of 

^r^ several tens of nanometers, a promising prerequisite for quantum sensing applications. 
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1. Introduction 

Tremendous progress in understanding and manipulating the nitrogen vacancy (NV) 
center in diamond throughout the last decade revealed its promising capabilities for 
quantum information and sensing applications. The basic foundations of coherent 
manipulation, high fidelity polarization and optical readout [1] paved the way for 
using the NV center as a fully controllable quantum bit capable of operating at room 
temperature with extraordinary long coherence times that may reach the millisecond 
range [2] . Whereas the coherent coupling and entanglement to nuclear spins of 
nitrogen [3, 4] and carbon- 13 [5, 6] has been demonstrated in numerous experiments, 
bringing diflFerent nitrogen vacancy centers to interaction remains challenging and has 
been demonstrated only recently [7, 8] . One approach for coupling distinct NV-centers 
makes use of their dipolar interactions [7, 9, 10], which is limited by the strong distance 
dependence of the coupling and therefore has been demonstrated only for very closely 
separated pairs. Another method consists of interconnecting the NV center solid state 
spin qubits with photons [8] , that has lead to extensive research in the design of cavities 
and photon couplings [11, 12, 13, 14]. In contrast to that, the coupling to phonons 
is much less studied. Whereas this mechanism serves as the prominent data bus for 
conditional quantum operations in the trapped ion approach to quantum computing [15], 
intrinsic phonon coupling is assumed to be inaccessible in macroscopic diamonds at 
room temperature. However, the coupling to magnetized nanomechanical oscillators as 
AFM cantilevers was successfully performed, allowing for the sensing of the vibrational 
mode [16, 17] and even for the coherent manipulation of the NV center electron spin 
state [18], that might provide the basic ingredient for future phonon mediated quantum 
networks [19]. Here we will show that a significant intrinsic phonon coupling can be 
expected in nanodiamonds at low temperatures. Those can be fabricated down to 4nm 
in size, additionally being capable of hosting NV centers [20]. We study the coupling 
strength to long wavelength (low frequency) acoustical modes and analyze the possibility 
to exploit these global modes for entanglement operations by creating a Raman-induced 
phonon coupling within the ground state electron triplet states of the NV center. 

2. Coupling of the nitrogen vacancy center to phonons 

The localization of the NV-center electronic states to the vacancy defect itself permits 
their description as superpositions of molecular orbit als, each associated with a dangling 
bond orbital of the atoms involved in the defect center [21, 22]. That is, according to 
the Csy symmetry of the defect center, four electronic states A4q\ = {|ai), |a2), le^), \^y)} 
can be constructed out of the dangling bond orbitals A^db = (l^^c)' Wc)^ Wc)^ Wn)} by 
linear combinations, with the index C referring to carbon and A^ to the nitrogen related 
bonds (see figure 1 (a)). 

In this framework the well known ground and excited state energy level structure 
of the NV~- center follows by associating the six electrons involved with the available 
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Figure 1. Energy level structure and phonon coupling (a) NV center energy 
levels and composition out of dangling bond orbitals. (i) NV center and dangling bond 
orbitals with symmetry Csy. (ii) Symmetric linear combinations of the dangling bond 
orbitals (LCAO) lead to the electronic states |ai), |a2), le^^), \ey), occupied with the six 
electrons (or equivalently two holes) of the NV center, (iii) Combining the electronic 
and spin wavefunction leads to the NV center energy level structure with spin singlets 
and triplets. Ground states correspond to two holes in the le^/y) state as illustrated in 
(ii) and excited states to one hole in |a2) and le^./^). Strain and therefore phonons cause 
displacements of the dangling orbitals (i), and consequently modify the energy levels in 
(ii) and (iii). The influence of the strain Sx = —(^xx on the excited state fine structure 
is illustrated in the inset, (b) Phonon coupling coefficient 77 (6) and frequency i^ of the 
lowest acoustical mode vs size (diameter) of a spherical nanodiamond. 



electronic states Aigi, taking the additional spin properties into account. Equivalently, 
and more simply, this can be described by two holes relative to a complete filling [21]. 
Calculating the coupling to phonons can be decomposed into two steps by noting that 
the long wavelength acoustical phonons considered here will introduce a periodic strain 
to the NV center. Thus in a first step we will discuss the effect of strain to the energy 
level structure and in a second step link the vibrational phonon mode to the strain 
property. For the former case we will follow the discussion presented in [21] (see also a 
similar discussion in [22]). The electron- nuclei Coulomb coupling can be described by 
an interaction of the form 



V = ^gi \cri){ai\ + ^hij{\fij\) \(Ji){(Jj\ + h.c. 



(1) 



i>j 



with i and j describing a summation over all possible A^db configurations and fij = fi—fj 
with fi/j the corresponding atom position vectors. Noting that the coupling coefficient 
depends on the relative distance of the atoms involved, it is obvious that stress related 
displacements Ui will influence the energy level structure and for small displacements, 
as expected by phonon effects, the strain perturbation on the level of dangling bonds 
can be described by 



5V 



E 



^hij{\fij\) \(7i){aj\ +h.c. 



(2) 
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with 5hij = V^hij 5u = V^hij ef the first order correction foUowing from an expansion 
of hij. Herein the tensor e^^ = du^/{dry) is related to the strain tensor e by 
ti^i, = l/2(e^^ + e^^). Considering a specific ground state to excited state triplet 
transition this results in the following Hamiltonian for the NV-center energy levels: 

//ei-phon ^ 25i \g){g\ + \e){e\{S, + S4 + O (3) 

wherein 5i = —({cxx + ^yy) and ^4 = (8/3'^ Czz with [21] ( c^ 610 THz and /3 accounting 
for the difference of the nitrogen coupling compared to the carbon related ones. I^') 
can be any of the triplet ground states, that exhibit an equal shift under the infiuence 



of strain, whereas for the excited states we find |e) ^ = ±C w4e^^ + [cxx — ^yyY ^r 
the states \Ex) and \Ey)^ respectively, and C = otherwise. For more details on the 
derivation of (3) we refer to Appendix A. 

Vibrational modes will introduce a time-dependent periodic strain to the system, 
that exhibits a simple description for long wavelength acoustical modes with wavevectors 
near the Brillouin-zone. In that case the discrete Bloch-type wavefunction, describing 
the displacement of the crystal lattice positions, can be described by a continuous 
displacement field in space u{r) [23] replacing the discrete lattice positions nhy r^ ^ r 
and the displacement for a specific mode a is thus obtained by [24] 



with k the wavevector, e^^^ the mode eigenvector, v{k) the angular frequency and a^ a^ 
the phonon creators and annihilators of the mode, respectively, and M the mass of 
the nanodiamond. For the linear acoustical modes near the zone center z/(/c) = c - k 
with c = 1.2 • lO^m/5 the speed of sound in diamond and k = 2 7t/1 out of applying 
periodic boundary conditions, with / the diamond length in the corresponding mode 
direction. Combining equations (3) and (4) allows for the calculation of the energy shift 
associated with the phonon coupling and exhibits the typical form of a deformation 
potential coupling describing local lattice compression and dilation as expected for the 
long wavelength acoustical phonon case. As an illustration let us explicitly give the 
electron-phonon coupling energy shift Hamiltonian for the specific case of a mode a with 
e^^^ II k pointing in x-direction and choosing \e) = IA2) as the excited state, favourable 
in the sense that it is not coupled to the singlet state via spin-orbit coupling [21, 25]. In 
that case the electron-phonon coupling takes the simple form 

^ei-phon = -r]v{-i) {a^ - a) |e)(e| (5) 

with the coupling coefficient 



that takes the role of the well-known Lamb-Dicke parameter. 

Adding a laser coupling with Rabi frequency ^ and frequency uol to couple one of 

the ground states to the excited state |e) and applying the canonical Schrieffer- Wolff 
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transformation U = exp{ir] {a'^ + a)|e)(e|) beside restricting the discussion to a single 
mode, leads to the following Hamiltonian in the rotating wave approximation [26, 27] (a 
discussion on how to include the mode relaxation and the finite excited state lifetime 
can be found in Appendix E) 



UJQ 



IT ^U , t , 

H = -— az + jya^a + 



-|e)(^|e-^^^^e^^(^^+^)+h.c. 



(7) 



2 

Herein a^ refers to the Pauli operator of the corresponding transition and ujq = cjo + t]'^ u 
with ujq the transition frequency and the second part following from the SchrieflFer- 
Wolff transformation. Note also that the phonon coupling mechanism is assumed 
to be completely originated from the electron-phonon coupling Hamiltonian of the 
corresponding transition, whereas the photon recoil based coupling mechanism has been 
neglected as it is orders of magnitude lower in macroscopic systems [26] . The coupling 
parameter rj for diflFerent nanodiamond sizes is shown in figure 1 (b) taking typical values 
of the order of 10~^ — 10~^ and decreasing with the nanodiamond radius oc r~^ for a 
spherical diamond of radius r. For general diamond shapes rj oc \/ljV with / the length 
in the mode direction and V the diamond volume, leading e.g. to a ry oc l/y^ scaling 
for a two-dimensional structure of area oc r^. The mode frequencies v for the lowest 
energy mode, for a spherical diamond also representing the frequency difference among 
neighbouring modes, is shown in figure 1 (b) taking typical values in the THz-range. 
These high mode frequencies are advantageous in the sense that neighbouring modes are 
well separated and additionally thermal occupation probabilities are low. The decreasing 
magnitude of both the coupling strength and the mode frequency limits the diamond 
size to several tens of nanometers. Hamiltonian (7) can be expanded in orders of the 
small coupling parameter r] with the zeroth phonon independent order providing the 
carrier and the first order the blue (oc a^\e){^g\) and red sideband (oc a\e){^g\) transitions. 
These phonon sidebands should be experimentally observable in low temperature 
emission spectra [28] ; however no such experimental investigations with nanodiamonds 
< 5 nvn have been performed to date. 

3. Phonon mediated gate interaction 

The phonon dependent sideband transitions offer the possibility of correlating the NV 
center state with the global vibrational one. This can be used as a source for gate 
operations [29, 30, 31] between different NV centers within the same nanodiamond. 
In one of its simplest versions a M0lmer-S0rensen type gate [29, 32] can be directly 
implemented on the ground-excited state transition, that has the advantage of being 
insensitive to the actual phonon state, therefore allowing gate operations even for 
thermal initial phonon states and circumventing the need for vibrational mode cooling. 
However such a concept suffers from the relatively short lifetime of the excited state 
limited by both optical decay to the ground state as well as non-radiative spin-orbit 
relaxations to singlet states. Moreover, in order to avoid off-resonant excitations to 
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Figure 2. Phonon mediated gate (a) Setup for creating a Raman induced phonon 
coupling within the ground state manifold by off-resonantly exciting the excited state 
with a carrier and blue-sideband transition. Dashed lines represent couplings that just 
appear in the double-path setup. Green lines are compensation pulses that become 
relevant for gates in the non-perturbative limit hi2 = ^/Ae < 1. (b) Final effective 
two-qubit gate interaction for the single (only continuous lines) and both-path (all 
interactions) configuration, (c) Effective gate frequency vs size (=diameter) of a 
spherical nanodiamond and different values of tvi and f<i2^ taking the coupling to the 
second neighbouring excited state into account. Herein Q2 = ^1^, ^1 = r]Q2 and 
hZ2 defined as in the text. Continuous, dashed and dashed dotted lines correspond to 
hZ2 = 0.1,0.05,0.35, whereas blue, red and green corresponds to hzi = 0.01,0.05,0.1, 
respectively. Horizontal lines (with triangle endpoints) denote the corresponding Feff 
value that imposes a limit on the maximal nanodiamond size shown for K:2 — 0.05 and 
K.2 = 0.1 by the vertical yellow lines. The area for which the ratio legate > Teff 
is coloured with decreasing saturation for increasing /^2- Note that this ratio is 
independent of the choice of ki for a given nanodiamond size. 



several excited states, the driving field strength has to be limited to values smaller than 
the typical energy gaps ^ AGHz^ therefore limiting the maximal gate speed. 

Here we provide an alternative concept that allows to perform the gate within the 
triplet ground state manifold. A Raman transition via the excited state provides the 
required phonon coupling, that otherwise would not exist according to the absence of a 
difference in the electron-phonon coupling among the ground state triplet states in (3). 
For this setup to work, a A-transition between ground and excited states is required, that 
for the NV center exists between the Is'+i), \g-i) electron-spin triplet ground states and 
either of the excited states \Ai) ^\A2) ^\Ei) ^\E2) ^ the latter forming equal superpositions 
of the rris = ±1 electron spin projections [21]. These A-transitions have already been 
successfully implemented and analyzed in experiments [33, 25], being accessible in the 
low temperature (< 10 X) and strain limit. Using circular polarized light allows to 
drive spin-selective transitions between those states that might be advantageous for 
tuning the detunings and couplings. Additionally the Raman-transition within the A- 
scheme has to be carefully tuned to a single sideband, that is either the red or blue 
one, as otherwise the phonon dependence of the effective ground state transition is 
cancelled by interference of those two paths. For this task the high phonon frequencies 
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of nanodiamonds are advantageous as they allow the single sideband addressing without 
stringent conditions on the Rabi- frequency of the coupling field. The gate interaction 
itself follows from a two-step process: First, the Raman transition provides a phonon 
dependent coupling within the ground state manifold. Second, this phonon dependent 
coupling can be used to implement a gate between two NV-centers by off-resonantly 
exciting the phonon state similar to the direct M0lmer-S0rensen approach. That way 
the excited state relaxation is suppressed by the off-resonance of the Raman transition, 
leading to an improved ratio between gate and relaxation time compared to a direct 
gate implementation on the ground-excited state transition as will be discussed below. 
Moreover, contributions of different excited states simply add up, therefore allowing for 
a straightforward integration of this effect into the formalism and not leading to an 
excitation of several states as might happen in the direct implementation. 

3.1. Setup and first effective form 

The setup for this gate is illustrated in figure 2 (a) where the dashed transitions are 
present only in the double-path (dp) setup, that leads to a complete ax (8) a^^-type 
coupling in the ground state manifold, whereas they are absent in the single-path (sp) 
one resulting in a gate in the reduced manifold Aii = {\g^i^g_i)^\g_i^g^i)}. The 
corresponding Hamiltonian in a frame rotating with the laser frequency for NV center 
k can be written as 

^^ = Y|e)(^+l|e-^* + ^|e)(^_l|e-(-+-)* + ^r7,^at|e)(^_l|e--* + h.c. (8) 

with Qi 0^ r]k^2- Herein the first and last contribution form the phonon-dependent 
Raman transition consisting of a carrier and blue sideband excitation, respectively, 
provided that ei ^ Qi and 62 ^ ^2 (the same can be achieved with a red sideband 
interaction as well). The second contribution describes the unavoidable carrier excitation 
associated with the sideband term and is not required for the gate interaction itself, 
however cannot be neglected either. Higher order terms in rjk have been omitted. 
For the double-path setup Hj^^ = H^ + ^^^L ^)^\q_i) with the second contribution 
corresponding to the first one by replacing \g+i) by \g-i) and vice versa, i.e. both of 
the couplings Qi and ^2 are present simultaneously on both transitions. Describing 
the off-resonance ratio between driving fields and detuning by /^i <C 1 and noting that 
e/e <C ^, the optimal choice of parameters is given by 1^2 "^ ^1 ^ and ci c^ 62 = l//^i ry f^2- 
Eliminating the off-resonant excited state results in the effective Hamiltonians 



^i,dp _ ^ , 

^eflf,/e — 2 ^ ^ 



ifsp 



. ^k 

'Y 



e^'^'*a^((T^ + l) + h.c. 
e''^^*aV+ + h.c. 



(9 a) 



(96) 



with Qk = i r2i (?7fc ^2) (ei + e2)/(eie2) and 5^p = 5^=2 and S'p = l/2S^=i with 
Sf, = 1/2 {Qj/ei + {-1T4 (1 + n) J^i/e2 + (-1)'^ Qi/(z/ + €3)). Moreover Ae = ei - 
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^2 + XVk^l/^^ where the last term accounts for the couphng induced shift of the mode 
frequency with x = 1/4 for the double and x = 1/8 fo^ the single-path configuration. 
The Pauli operators a^^, a^ and a+ are defined in the ground state manifold (Is'+i), \g-i)} 
and n denotes the phonon number operator. The terms involving phonon excitations 
correspond to paths involving both the sideband and carrier transition whereas non- 
combined paths lead to single transitions or AC Stark shifts not associated with an 
eflFective phonon (de-) excitation, respectively. Note that the last contribution in 5^ 
originates from the Q2 carrier interaction and can be significantly larger by a factor of 
I/t] than the preceding terms. Interestingly, taking into account the coupling to several 
excited states, contributions arising from \Ai) and IA2) (as well as \Ei) and |£'2)), both 
equal m^ = +1 and m^ = — 1 superpositions, have opposite signs and therefore subtract. 
This limits the maximal amplitude of the uncorrelated flip contributions arising from 
carrier transitions to hif A with A ^ 4 GHz the excited states energy splitting, that can 
be considerably lower than expected from the coupling to a single excited state. 

3.2. Second effective form and gate Hamiltonian 

In a second stage we consider Hamiltonian (9a, 96) for two NV centers (k=l,2) and 
choose Ae = l//^2 ^ with /^2 <^ 1 denoting the off-resonance of the corresponding 
transition. That way the phonon transition is only virtually excited which allows to 
obtain a second effective form 



k=l,2 
( 

'2 "^ 4 2 



esp -j Q 

k=l,2 

corresponding to the flnal gate Hamiltonian. Herein the effective gate frequency is 
deflned as Q^^te = ^1 ^2/Ae, ^f = (^f + (f^i + f^2)f^fc/Ae and ^^^ = 5^^ + ^^^g^. 
For the double-path scheme this corresponds to a a^; (8) a^^-type gate rotating the 
states within the two-qubit manifolds A^i = {l^'+i, ^'-i), l^'-i, ^'+1)} and A42 = 
{\g^i^ g^i)^ \g_i^g_i)} at a Rabi-frequency f^gate ^ ^1 ^2 (^^2)- Additionally there exist 
uncorrelated single qubit flips between the ground state levels with 5^^ ^ /^i r]Q2 that are 
by an order l//^2 larger than the gate term itself (see flgure2 (b)). However both terms 
commute what allows to remove the single qubit flips by a simultaneous echo-vr-pulse 
in az or ay on both NV centers, leaving the gate interaction unchanged but adding a 
negative sign to the uncorrelated single-flip contributions. Interestingly, choosing ay for 
that task offers the possibility to decouple the ground state system from decoherence 
processes as well, therefore extending the coherence time signiflcantly. To conclude, a 
pure gate interaction can be achieved by adding any periodic pulsed decoupling sequence 
in ay acting on an inter-pulse timescale larger than the one required for the effective 
Hamiltonian form to be valid, i.e. At ^ {^^)~^ ^ (^^2)""^- For a two-qubit 7r/2- 
rotation this gate interaction including the echo-refocusing is illustrated in flgureS (a). 
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The single-path scheme behaves in a similar way, with the gate interaction restricted 
to a rotation in the Aii manifold, more challenging with respect to the initial state 
initialization on the nanometer scale of adjacent NV centers that requires individual 
addressing. For the case of equal a^-contributions, i.e. identical configurations on both 
of the NV centers, the gate and single qubit contributions commute again and the latter 
can be removed by an echo pulse and combined with decoupling sequences in ax and 
cr^, analogue to the two-path situation. Uncorrelated transitions do not occur as they 
do not form sideband independent paths. Note that dependent on the parameters, the 
AC Stark shift contributions might be an obstacle in adjusting the off-resonance for the 
second effective gate Hamiltonian form. Therefore, choosing smaller Rabi frequencies 
{Q2 ^ V ^^) might be advantageous to suppress the predominant infiuence of the carrier 
term. 

3.3. Comparison to a direct gate implementation 

At this point it is interesting to compare the relevant timescales of this Raman-induced 
scheme to a direct gate implementation on the ground-excited state transition. For 
the direct gate implementation the conditional gate operation is directly implemented 
between one of the ground and the dissipative excited state, i.e. the Raman transition 
step is omitted. In the Raman-induced implementation the gate frequency oc /^i /^2(^^2) 
and has to be compared to the effective decay rate Feff = /^^F with F c^ 15 MHz the 
excited state decay rate [33] suppressed by the probability /^^ of actually populating the 
excited level. Q2 is limited by the off-resonance to the carrier transition as described 
above. The corresponding ratio follows as r^gate/rgflp ^ {1^2/1^1) (^^2/r). In contrast 
to that a direct gate implementation leads to a gate frequency f^g^f oc tviijQ with 
the same Q limitations, that in this setup has to be related to the bare decay rate F, 
therefore f^g^f /F = /^i rjQ/T. Note that this is by a factor /^ <C 1 worse than for the 
Raman-induced gate scheme. 

3.4' Time- conditioned gate 

To improve this ratio larger /^2 values are advantageous, corresponding to smaller off- 
resonances with respect to the intermediate sideband states. Interestingly the condition 
/^2 <C 1 can be significantly relaxed for the double-path setup, noting that in this 
case the time evolution following from Hamiltonians of the form (9 a) can be exactly 
integrated [32] as will be shown in Appendix B. An important prerequisite in that 
regime consists of compensating the phonon number dependent terms appearing in (9 a) 
(the 'ry^-terms') to ensure the commutativity of the gate relevant term (the (ik-^eim) 
with the contributions that do not involve an effective phonon (de-) excitation (the 
5j[;^-contribution) . Such a compensation can be achieved by the green compensation 
couplings illustrated in figure 2 (a) and also leads to significant improvements for the 
gate in the perturbative regime (10a), (106) in cases when it is implemented not deep 
within the K.2 ^ ^ limit. Replacing the condition Ae <C 1 by the gate time 
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condition tgate = m2Ti//\e (m G IN), the resulting evolution can still be described 
by Hamiltonian (10a), even this does not hold for intermediate time-steps. That way 
K2 = \JQ ji^T^rn) with Q the gate rotation angle taking the value K2 = l/(2v/2) for 
creating a maximally entangled state {9 = 7r/2) with tti = 2, providing that the phonon 
population refocuses before the intermediate echo pulse is applied (see figure 3 (b)). 
A similar gate that allows one to perform the gate in the non-perturbative regime despite 
maintaining its independence on the phonon state can be constructed out of the single- 
path configuration by adding a continuous microwave driving within the ground state 
triplet manifold [30] and we will discuss that idea in Appendix C. 
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Figure 3. Simulated gate interaction Simulation of the double-path gate 

interaction (10a) for a 15 nm diamond and (a) ni = 0.05, i<i2 = 0.05 and (b) /^i = 0.05, 
K,2 = 1/(2 \/2) — 0.35 (time-conditioned version). The upper figures show the 
population of the states |^+i,^+i) (blue), \g-i,g-i) (red) and |^-i,^+i), |^+i,^-i) 
(black). Dashed lines illustrate the interaction including the uncorrelated flips that are 
refocused by an echo pulse, whereas continuous lines focus on the pure gate interaction. 
The lower plots show the phonon population during the gate interaction that can be 
significant for the time conditioned (close-resonant) gate. 



3.5. Size limitation 

Figure 2 (c) compares the gate speed to the excited state decay limitation for different 
configurations including the time conditioned 1^2 = 0.35 case. The maximal size 
limitation is determined by the ofl[-resonance parameter /^2 and is given by < 25 nm 
for the /^2 = 0.05 and < 35 nm for /^2 = 0.1. Here it should be noted that the lifetime 
in nanocrystals is increased compared to the bulk counterpart due to the strong change 
in the refractive index and dielectric screening effects [20, 34] with typical values [35, 36] 
Fnd — 1/2 r = 7.5 MHz. Therefore longer excited state lifetimes have to be expected 
what was not taken into account in figure 2 (c) as the exact magnitude is very sensitive 
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to the substrate environment [34] . Provided that the Rabi frequencies can take the same 
values as in the bulk counterpart, a decay rate halved in value would double the ratio 
Qgate/^effi therefore making the conditional gate operator less prone to spontaneous 
decay for a given diamond size. For a spherical diamond this would increase the maximal 
radius by a factor of ^/2. 

3. 6. Influence of dipolar couplings 

Up to now we neglected the influence of dipolar couplings, that is the optical dipolar 
coupling [37] on the ground-excited state transition as well as the magnetic equivalent [9, 
10] within the ground state manifold. Whereas this is a good approximation for speciflc 
conflgurations, e.g. an angle of 54.7° of the axis connecting two NV centers to an 
equally oriented symmetry axis for which the dipolar couplings are exactly zero, for 
other conflgurations they can be of the same magnitude as the gate interaction itself. 
These couplings take values of jopt — Stt • 52AMHz (optical ground excited state 
coupling) contributing as jopt ^ ^liopt due to the excited state off-resonance, and 
imag — 27r • lOA kHz (magnetic ground state coupling) for an NV center distance of 
r = 10 nm. That is the effect of the dipolar interactions cannot be neglected in 
the general case. A detailed discussion about how to include the effect of dipolar 
couplings in the formalism is presented in Appendix D. As a result this requires to 
replace the detunings by Ck ^ ek — iopt/2 and the gate frequency for the both-path 
setup by f^g^, -^ f^g^, - 2 jopt - imag/2 and Q^l^^ -^ Q^l^^ + imag/2 for the Mi- 
and A^2-interaction, respectively, as well as f^gate ^ ^gate — 4 jopt for the single- 
path setup. Herein jopt ^ ^liopt- Therefore, as long as jopt <C ly together with 
identical coupling conflgurations on both NV centers, both coupling mechanisms can 
be combined by taking into account the modifled detuning conflgurations and adjusting 
e/e correspondingly. 

4. Summary 

In summary we analyzed the coupling of nitrogen vacancy centers to long wavelength 
acoustical phonons, a mechanism capable of mediating gate interactions between NV 
centers for nanodiamond sizes of several tens of nanometers. Exploiting the existence of 
a A-scheme in the low temperature and strain limit, fully noise decoupled two qubit gates 
can be constructed within the ground state manifold, even in the presence of dipolar 
couplings. This might be interesting for the creation of entangled states but also for 
manipulating the phonon mode itself, that is the control of the motional degrees of 
freedom, e.g. the cooling of vibrational modes. Moreover, the realization of entangled 
states in nanodiamonds could have crucial application for future sensing protocols. 

Author ^s note: While flnalizing this manuscript we became aware of a 
similar investigation [38] studying phonon induced spin-spin interactions in diamond 
nanobeams. 
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Appendix A. Calculation of the electron-phonon coupling [21, 22] 

The effect of strain induced by phonons, i.e. the electron-phonon coupling, can be 
calculated by analyzing the effect of a change in the intra-atom distance on the Coulomb 
coupling interaction equation (1). Assuming that this displacement is small, a realistic 
assumption by restricting the analysis to the long wavelength acoustical modes or more 
precise to the one defined hy kl = n2 7r with n = 1, for which the wavelength is 
given by the length scale of the diamond crystal, it suffices to expand the coefficients 
hij to first order in the atom displacements Ui out of the equilibrium positions fj^ (i.e. 
f^^fi = f^ + Ui). With the additional spherical symmetry assumption that the energy 
will just depend on the absolute value of the relative two atom displacement, that is 
hij = hij{\fi — rj\) one obtains hij = h^j + 5 hij with 

*« = i5ir^^W^ (f;°-^°) («.-«,.) (A.i) 



ri ^j 



1 



dhij 


{n 


-rj) 


d 


n- 




d 


Ti - 


-rj\ 



\' i ' j \ ^ \' I ' j\ 



(rt - n°) e {rf - r!) (A.2) 



where in the second line we introduced the displacement tensor e^^ = dUfj^/drj, and the 
derivation is evaluated at the equilibrium position of the atoms. Inserting the explicit 
expressions for the positions r^ as defined in figure 1 (a) and using the expressions 
of the electronic states A^ei in terms of the dangling bonds orbitals A^db, the strain 
Hamiltonian (2) can be rewritten as 

^Vm^x = -2Ce^^|e^)(e^| - 2Qeyy\ey){ey\ - 8/3^ C^^^ l<^2)(<^2| . . 

- C {^xy + e^^) (|e^)(e^| + h.c.) . 



Herein we defined 

c = 



2 dhij{\ri-rj\) 



3 d \fi — fj\ 



q (A.4) 

o,c 

with the index C referring to the coupling for two carbon atoms and the difference for 

the carbon-nitrogen case is accounted for by the factor /3. The quantity q denotes the 

next neighbour distance in the diamond lattice and is equal to g = ^3/8 {r^ — fj). 

Moreover couplings between \a2) and \ex^y) levels have been neglected, justified by the 

large energy separation. Those would correspond to strain induced transitions between 

the ground and excited states of the NV center and consequently do not play a significant 

role. Note also that the energy level \ai) has been neglected as it is delocalized in the 
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valence band and does not contribute to the properties of the NV center energy level 
structure. 

Out of equation (A. 3) the impact of strain on the NV center energy levels can be 
calculated by using the expression of the energy levels in terms of the electronic states 
Aie\ as provided e.g. in [21]. For the 'two hole' description the strain perturbation 
Hamiltonian takes the form 



H^ 



el— phon 



m 



Me 



lei + lei ® SVj 



Me 



1 



spin 



(A.5) 



with lei the identity on Aiei and Igpin the one on the spin degrees of freedom. Projecting 
Hamiltonian (A.5) on the NV center energy level states finally leads to 



^eT-phon = 2^1 {\go){go\ + \9+l){9+l\ 4 

for the ground state electron spin triplet states and 

f S1 + S4 
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(A.6) 
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(A.7) 



/ 



for the excited state triplet states in the basis {|^i), 1^2), \Ex)^ \Ey)^ l^i)? 1^2)} with 



^1 = -({^xx + eyy), S2 = -({e^x - e. 



yyj'> 



s. 



S \^xy ~r ^yx) 



and S4 = — 8/3^C^2 



Combining the Hamiltonians (A.6) and (A.7) and neglecting the off-resonant couplings 
between {|^i), 1^2)} ^ (l^i)? 1^2)} results in Hamiltonian (3) of the main text. 

In the case of phonons the displacement can be expressed in terms of Bloch type 
wavefunctions, such that for a specific mode a, and noting that in the framework 
of classical elasticity theory applicable in the long wavelength limit the microscopic 
structure can be replaced by a continuous displacement field (r^ -^ r) , the displacement 
takes the form 



Un{r) 



h 



2Mu{k) 



e (-) 



Akf 



a«e'"" +a'e 



i kr 



(A.8) 



with M the total mass of the system, e ^^^ the eigenvector of mode a, k the wavevector 
and iy{k) the angular frequency and a^ aj^ the mode annihilator and creator operators, 
respectively. With the definition of e^^ = dui^/dr^y and using that fcr <C 1 this 
results in equation (4) in the main text, that together with the electron-phonon coupling 
Hamiltonians (A.6) and {A.7) complete the analysis of the phonon infiuence on the NV 
energy level structure. 



Appendix B. Exact integration of the gate Hamiltonian 

The time evolution following out of Hamiltonian (9 a) can be integrated exactly for 
commuting state operators (e.g. for the double-path gate setup), what allows to 
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overcome the ^2 <C 1 limit by identifying appropriate time conditions. We will consider 
the generic form 

H{t)=i{^{t)da^ -Y{t)d^a) (B.l) 

with [O, O^] = 0. For the double-path setup considered here (9a) 



7W = fe^^'^* 



and 



(B.2) 
(B.3) 



Using the properties of the displacement operator D(a) = exp(«a'^ — a*a), the time 
evolution of (B.l) follows as | 



with 



[/ = re-^/^(*')^*' = D(«(t)C')exp(^ \^(3{t)6& -/3*{t)d^d\] (B.4) 



a{t) 



m 



/dt'7(0^4|(e--l) 



*dt'7(0 r dt"^*it")=t^^ it + ^ [e^^^^-l] 



(B.5) 



(B.6) 



/o JO 4Ae V Ae 

Noting that the phonon dependence only appears in the displacement operator, it can 
be eliminated by choosing the gate time tgate such that Aetgate = m ■ 2% with m E 1j 
in which case tt(tgate) = 0. Implying that condition is fulfilled, the total evolution 
corresponds exactly to the one out of (10 a) (up to local contributions of the first order 
effective Hamiltonian and global phases), i.e. 



U{tgate) = exp 



n' 



2Ae 



ialal + 2 [al + al]) 



t 



gate 



(B.7) 



Appendix C. Microwave assisted gate 

A gate interaction with similar properties as the double-path gate discussed in the 
main text can be constructed out of the single-path setup combined with a continuous 
microwave driving of the states l^'o) ^ Is'+i) ^^d l^'o) ^ \9-i) (see figure CI (a)). Mainly 
this allows to perform the gate in the non-perturbative regime implying a time condition 
to assure the independence of the phonon state. As a side eflFect such a driving also 
decouples the gate from ground state decoherence. 

Considering the single path setup alone, the gate term in {9b) includes the operators 



cr± = - (cr^ ±icry) • 



(C.l) 



In that case a closed phase space trajectory, that is a refocusing to the initial phonon 
state at a specific time, is prevented by the rotation around two orthogonal axes {ax 
and ay) [30]. In a more formal way the non-commutativity of ax and ay prevents the 
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Figure CI. Microwave assisted gate: (a) Gate setup. The gate consists 
of a single blue sideband transition connecting \g-i) to l^'+i) (realized by a Raman 
transition via the excited state) and a continuous microwave coupling fluw between 
the states |^o) and |^±i). An additional compensation coupling (green line) removes 
the phonon number dependent ac- Stark-shift, (b) Simulated gate process for CImw = 
10 MHz, Qi = 7/^2 = Q.3MHz, ei = OAGHz, Ae = b3MHz and a diamond size 
of 15 nm. The time condition is given by tgate = 47r/Ae (n=2). For simplicity and 
clarity the carrier contributions of the blue sideband laser have been neglected in the 
simulation. The cyan dashed line corresponds to the sum of all populations including 
the state |^o)- (c) Possible implementation of the gate process. Blue lines denote a 
time evolution under the Hamiltonian form (C. 2) and ^e G {ctP"^, 5'^, 5'^}. 



exact integration of the gate Hamiltonian as described in Appendix B. Now adding 
a continuous microwave driving such that the ay contribution is suppressed, removes 
those difficulties, and additionaUy leads to a {ax cra;)-type gate similar to the double- 
path gate proposal that rotates states both within Aii and A^2- One possibility to 
achieve that task would be to continuously drive the states \g^i) ^ Is'-i) [30, 31], 
what however is not really practicable for our setup. Therefore we will incorporate 
the full ground state triplet and show that a driving of the form l^'o) ^ \9±i) will 
be suitable for this task as well. In the following we will refer to the qubit operators 
within (Is'+i), |fl'-i)} by the indices 'pm' and denote the Pauli spin-1/2 operators in that 
manifold as {(jP"^, a^"^, a^"^, 1pm} whereas we will denote the spin-1 operators in the 
ground-state triplet manifold as {Sx^ Sy^ S^^t}. 

The Hamiltonian of the total system including the microwave driving Hmw is given 
by 



H 



EztO) , ztO) , ztO) , o-O) 
^MW ^ -"ipm ^ -"cr?'" ^ -"gate 



(C.2) 



=1,2 
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Q=-Qi {r]Q2) 
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where we assumed that the ry^-terms have been successfuUy compensated to 
avoid the existence of phonon-number dependent terms (see green laser couphng in 
figured (a)). This describes the microwave driving, the global ac-Stark shift of the 
\g±i) states with respect to the l^'o) state, the relative shift of the \g±i) states and the 
gate relevant term, respectively. Note that the global H^^^ term can be neglected if the 
microwave frequency is already tuned to the ac-Stark shifted states. 

In an interaction picture with respect to the continuous microwave driving and 
assuming that ^mw ^ {^i^sp^^i}^ the following substitutions can be made by 
neglecting fast rotating terms ('rotating-wave approximation') 



(C.4) 



^x = ^x^x ~ ^y^y 


^ 
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!, the Hamiltonian takes the form 
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(C.5) 



That way the a^-contribution is suppressed, all contributions commute, ideally the 
state \go) is never populated in the interaction frame and importantly the operators 
appearing in the gate Hamiltonian part (the second term in (C.5)) commute what allows 
to integrate the time evolution exactly as described in Appendix B (corresponding to 
O = 1/2(3/4(7^"^ — 1/4 1pm) in equation (B. 3)). Note that, as already discussed in 
the double-path setup, the uncorrelated single flip interactions can be removed by a 
single echo vr-pulse in a^"^ or a^"^. Therefore, including an echo-pulse the effective time 
evolution is exactly given by 

tgate I With il 



Uit, 



exp —I 



^gate pm 
2 ^ 



(J, 



pm 



gate 



if the time condition Ae t, 



gate 



choice with respect to the gate time n = 2 {n 



8 8 A. «=■«' 

n (27r) (n G IN) is fulflUed. That is, for the optimal 
1 cannot be realized due to the required 
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echo pulse) this leads to /^2 = ^/Ae = ^(4/3)^ O/tt for performing a two qubit rotation 
^gate • tgate = ^- Here it is interesting to note that the (maximal) ratio between the gate 
speed and the effective excited state decay rate is independent of the absolute values 
of the laser Rabi frequency (independent of /^i = ^2/^ = ^i/^i as defined in the main 
text) and therefore the condition ^mw ^ ^ does not alter the maximal nanodiamond 
size limitations as depicted in figure 2. However the absolute magnitude of the gate 
speed (^ /^i) decreases with a decreasing microwave field, such that other limitations 
as the T2-time of the ground state triplet states can replace the effective excited state 
decay rate as the limiting quantity. 

Note also that the continuous microwave driving decouples the gate interaction 
from ground state decoherence, an effect that can be modelled as a fiuctuating energy 
shift i/decoh = 5{t)/2Sz and is suppressed according to equation (C. 4), such that the 
limiting quantity will be given merely by Ti for a strong driving (for a more complete 
discussion of the decoupling method we refer to [9]). As a final remark, the interaction 
frame with respect to X^j^i 2 ^mV ^^^ ^e implemented using the pulse sequence (see 
figured (c)) 

exp {-I H,,, t) = U-.,,e-'''"^Ul,e-''''l' (C.7) 

with H given by(C.2) and Uint = exp^—in Se) and Se G {a^'^^ Sy^ Sz} (noting that 
UintSxU^^^ = —Sx whereas Hint is invariant under this pulse sequence). 

Appendix D. Influence of dipolar couplings on the gate interaction 

On the nanometer distance between NV centers, dipolar interactions can play a 
significant role and might itself provide the conditional coupling interaction at the same 
time disturbing the phonon induced mechanism. The optical dipolar interaction on the 
ground-excited state manifold can be described by [37, 39] 

^opt,dip = -^ (|e,^+i)(^+i,e| + |e,^_i)(^_i,e| + h.c.) (D.l) 

with the coupling constant [39, 40] (n /cq r <C 1) 



-^'"P* " 2 {nkoly ^^' • ^2 - 3 (pi . er) {h • e.)) 

~ 27r ■ 52.4 MHz f j (j5i ■ p2 - 3 {pi ■ e,) (p2 ■ e,)) 



(D.2) 



with r the distance between the NV centers, V the spontaneous decay rate (F = 
15 MHz [33, 41]), fco = 27r/Ao the vacuum wavevector of the transition and Aq = 637 nm^ 
n the refractive index (n = 2.4), ^^ the fraction of emissions into the zero phonon line 
(<^o — 0.03 [42]), Pi the normalized dipole moment direction of NV center i and e^ the 
unit vector in the direction of the axis connecting the two NV centers. Additionally 
there exists a magnetic dipolar interaction within the ground state triplet manifold, 
that, nevertheless being orders of magnitude weaker than the optical one, can play a 
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role for far detuned gates. It takes the form [43, 9] 
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(D.3) 



wherein /xq denotes the magnetic permeabihty, 7e^ the gyromagnetic ratio of the electron 
spin and a^ the Pauli z-matrix defined within (l^'+i), \g_i)}. 

For the derivation of the first effective Hamiltonian form (9 a, 96) the infiuence of 
the magnetic dipolar interaction on the excited state elimination can be neglected, and 
for the double path scheme the following Hamiltonian is obtained assuming identical 
configurations on both NV centers 
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This effective form is valid for any magnitude of the dipolar coupling as long as 
|r^i| c^ \r]Q2\ <^ |(iopt/2) — ei\ and additionally \i^l{j/^k)\ <^ 1 for \j\ > \ek\ to avoid 
a quasi resonant second order process to the double excited state \ee). Compared 
to the derivation without dipolar interaction (9 a), the dipolar modified Hamiltonian 
follows essentially by replacing the detuning Ck ^ ek — (j'opt/2) what has a clear 
interpretation as a coupling to the dressed states of the optical dipolar interaction (see 
figure Dl). Additionally direct off-resonant ly suppressed dipolar interaction terms of the 
order jopt ^ jopt ^1 appear in the Hamiltonian together with the magnetic equivalents. 
Interestingly, for the case of equal couplings on both NV centers, only couplings to the 
1+) dressed state occur (except for the state independent phonon terms described by 
Qa)i explaining why the detuning Ck + (iopt/2) is not a relevant quantity in the effective 
form (D.4). This is originated in the fact that |— ) couplings cancel by interference due to 
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the different sign of both paths as illustrated in figure Dl. However note that this is only 
valid for identical coupling configurations on both NV centers, otherwise replacements 
of the form (here for k=l, upper indices refer to k) 

^1 ^2 ^1 ^2 ^ I 1 iopt / ^1 ^2 I ^1 ^2 \ /T-v n\ 

( ■ /o\9 ,9 "^ O O I / • /o\9 ,9 "^ / • /o\9 ,9 I \^'^) 



are required in 5^^, ^/^ and f^^, taking into account couplings to both dressed states 
with the first part describing the dipolar independent and the second one the dipolar 
induced contributions. The terms in jopt always involve terms of both fc = 1 and A; = 2 
such that the replacements are obvious in that case. 

In the absence of the magnetic dipolar term in(D.4) all terms commute and 
therefore the optical dipolar contribution just adds to the phonon induced gate 
interaction that can be derived from the first two terms as in (10a). A small modification 
of this behaviour arises from the magnetic dipolar contributions, that do not commute 
with the (Ja^-terms in (D.4). However due to jmag <^ (^^^^ ^) the magnetic dipolar term 
contributes to (D.4) in the secular approximation as 

^mag,dip-^^(a>.^ + ajaj) (D.7) 

that commutes with all other contributions. Thus, the second effective Hamiltonian 
form analogue to (10 a) is given by 

\ J k ^^^ (D.8) 

I -^QP^ /^l ^2 I ^1 ^2 I ^1 ^2\ I Jmag / i 9 , ^1 ^2\ 
+ -Y \^^ ^^ + ^y^y + ^z ^z) + ^p y^z ^z + ^y ^y) 

wherein the single fiip-contributions (the first term) can again be removed by an echo 
pulse and all contributions commute. The a^a^-contributions just correspond to a global 
phase in the A4i and A^2-nianifolds such that in total the gate frequency is given by 
f^gate = -^V^^ + 2iopt + imag/2 for a Ml rotation, and legate = -^V^^ " imag/2 for 
a Ai2-rotation, respectively. 

A similar analysis can be carried out for the single-path setup leading to the same 
replacements of the detunings by e ^ e — (jopt/S) in the uncoupled form {9b) and direct 
dipolar contributions as in (D.4) with 



^ _ 1 

Jopt "^ Jopt 



ni r]'^ni{i + n) a 



In total the gate interaction arises from replacing f^gate by f^gate — 4j'opt in (106) beside 
the detuning replacements. 

Appendix E. Mode decay and excited state decay 

In here we will briefiy comment on how the spontaneous decay and the coupling among 
vibrational modes can be incorporated into the gate formalism. Both effects can be 
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Figure Dl. Dipolar coupling (single flip process): Illustration of a single 
flip process of NV center 2 from |^+i) to \g-i) in the presence of optical dipolar 
interactions. The initial and final states are coupled to the dressed states |+)(+i) and 
|— )(+i) here illustrated for the case |jopt| > kl (assuming ei c^ 62 — e) and e < (red 
detuned). Blue lines denote couplings that are independent of the dipolar interaction, 
whereas red dashed ones describe dipolarly induced processes. Upper indices denote the 
corresponding NV center and lower ones refer to a specific coupling type. Paths related 
to the |— )(+i) state cancel by interference due to the negative sign if 1]!^ ^ = fl^^\ i.e. 
in general if the configurations on both NV centers are identical. Analogue processes 
can be identified for other single fiip and |^+i,^_i) ^ |^_i,^+i) conditional coupling 
processes. 



described using a master equation approach. For the spontaneous excited state decay 
this takes the form 



dp 



(2a_p(7+- {p,a+a_}) 
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with cr+ = |e) {{g+i\ + {g-i\) and a_ = (cr+) , F the excited state decay rate, {a, 6} = 
ab + ba denoting the anti-commutator and a± = a±e^^^^^ '^^K The phonon mode 
relaxation can be modelled as [26] 
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with jy the mode frequency as defined in the main text, Q the corresponding mode 
quality factor and nth the thermal phonon population. Both of those equations hold on 
the level of the original Hamiltonian after the Schrieffer- Wolff transformation, i.e. in 
the same frame as Hamiltonian (8). 

The elimination of the excited state manifold, i.e. calculating an effective 
Hamiltonian as in (9 a, 96), corresponds to a unitary transformation T [44] such that 
Hl^ = VgT HT^Vg with Vg the projector on the ground state manifold and H{t) the 
original Hamiltonian as e.g. given by (8). Up to second order the transformation is given 
by [44, 45] (for the case of a Hamiltonian that is purely off-diagonal, e.g. only includes 
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terms that couple the ground to the excited state) 

T = l+iS with S = f H{t') dt' (E.3) 

with H{t') the ground-excited state couphng Hamiltonian [H^^ (8) for the single path 



jdp 



setup or H^^ for the double-path equivalent), and the effective Hamiltonian follows as 
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An arbitrary operator O in that effective frame can therefore be calculated by 
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which for the specific case of the ladder operators a_ and a results in 

eff,I 



gZ''' = —i CF- S, aeffj = a — - [S, [S^ a\] 



(E.4) 



(E.5) 



(E.6) 



eflf,I 



eff,I 



and a^'^ = fcr!_"'M and a^^j = (aeff,i) . Therefore the evolution in the first 
effective frame is correctly described by replacing the operators appearing in the master 
equations (E.l) and (E.2) with the effective ones as defined in (E.5, E.6). Additionally 
one has to account for the ground state dephasing by either a stochastic approach 
or a master equation approach in the Markovian limit [9] , unchanged by the unitary 
transformation for the first effective form, as T is diagonal by construction within the 
ground state manifold. 

To give an explicit example we will consider the double-path setup and assume 
equal coupling configurations rji = r]2 = ry, which results in (with a^;, 1 defined within 
the ground state manifold, e.g. ax = \g^i){g-i\ + h.c.) 



a5^ = (l + a,) fi^e— ^* + 



flo 



«eff,I = Ct 



2ei 



2(62 + U) 

2ei 



e-(e2+-^)* + ^at^e— 2* 



11^2 



r]^2 



a 



2e2 

1^2 



^ii^t 



(E.7) 
(E.8) 



262 " 2(e2 + z^) 

Out of (E.7) it is obvious that the decay rate F appearing in (E.l) can be replaced by 
an eflfective one of the order Tes,i ~ i^i as expected by the fact that this corresponds to 
the probabihty of populating the excited state. 



Coupling of nitrogen vacancy centers in nanodiamonds by means of phonons 22 

References 

[1] Jelezko, F., Gaebel, T., Popa, I., Gruber, A. & Wrachtrup, J. Observation of coherent oscillations 

in a single electron spin. Phys. Rev. Lett 92, 076401 (2004). 
[2] Balasubramanian, G. et al. Ultralong spin coherence time in isotopically engineered diamond. 

Nature Mater. 8, 383-387 (2009). 
[3] Gaebel, T. et al. Room-temperature coherent coupling of single spins in diamond. Nature Phys. 

2, 408-413 (2006). 
[4] van der Sar, T. et al. Decoherence-protected quantum gates for a hybrid solid-state spin register. 

Nature 484, 82-86 (2012). 
[5] Dutt, M. et al. Quantum register based on individual electronic and nuclear spin qubits in diamond. 

Science 316, 1312-1316 (2007). 
[6] Neumann, P. et al. Multipartite entanglement among single spins in diamond. Science 320, 

1326-1329 (2008). 
[7] Dolde, F. et al. Room-temperature entanglement between single defect spins in diamond. Nature 

Phys. (2013). 
[8] Bernien, H. et al. Heralded entanglement between solid-state qubits separated by 3 meters. 

arXiv:1212.6136 (2012). 
[9] Albrecht, A. et al. Self-assembling hybrid diamond-biological quantum devices. arXiv: 1301. 1871 

(2013). 
10] Bermudez, A., Jelezko, F., Plenio, M. B. & Retzker, A. Electron- mediated nuclear-spin interactions 

between distant nitrogen- vacancy centers. Phys. Rev. Lett. 107 (2011). 
11] Wolters, J. et al. Enhancement of the zero phonon line emission from a single nitrogen vacancy 

center in a nanodiamond via coupling to a photonic crystal cavity. Appl. Phys. Lett. 97, 141108- 

141108 (2010). 
12] Riedrich-Moller, J. et al. One-and two-dimensional photonic crystal microcavities in single crystal 

diamond. Nature Nanotech. 7, 69-74 (2011). 
13] Faraon, A., Santori, C., Huang, Z., Acosta, V. & Beausoleil, R. Coupling of nitrogen- vacancy 

centers to photonic crystal cavities in monocrystalline diamond. Phys. Rev. Lett. 109, 33604 

(2012). 
14] Hausmann, B. et al. Integrated diamond networks for quantum nanophotonics. Nano Lett. 12, 

1578-1582 (2012). 
15] Leibfried, D., Blatt, R., Monroe, C. & Wineland, D. Quantum dynamics of single trapped ions. 

Rev. Mod. Phys. 75, 281 (2003). 
16] Kolkowitz, S. et al. Coherent sensing of a mechanical resonator with a single-spin qubit. Science 

335, 1603-1606 (2012). 
17] Arcizet, O. et al. A single nitrogen- vacancy defect coupled to a nanomechanical oscillator. Nature 

Phys. 7, 879-883 (2011). 
18] Hong, S. et al. Coherent, mechanical control of a single electronic spin. Nano Lett. 12, 3920-3924 

(2012). 
19] Rabl, P. et al. A quantum spin transducer based on nanoelectromechanical resonator arrays. 

Nature Phys. 6, 602-608 (2010). 
20] Tisler, J. et al. Fluorescence and spin properties of defects in single digit nanodiamonds. ACS 

Nano 3, 1959-1965 (2009). 
21] Maze, J. et al. Properties of nitrogen-vacancy centers in diamond: the group theoretic approach. 

New J. Phys. 13, 025025 (2011). 
22] Doherty, M., Manson, N., Delaney, P. & Hollenberg, L. The negatively charged nitrogen- vacancy 

centre in diamond: the electronic solution. New J. Phys. 13, 025019 (2011). 
23] Ashcroft, N. & Mermin, N. Solid State Physics (Brooks/Coole Thomson Learning, 2006). 
24] Yu, P. & Cardona, M. Fundamentals of semiconductors: physics and materials properties (Springer 

Berlin, 1999). 



Coupling of nitrogen vacancy centers in nanodiamonds by means of phonons 23 

[25] 
[26] 



[27; 
[28; 
[29 
[30 

[31 
[32; 
[33; 
[34; 
[35; 

[36; 
[37; 
[38; 

[39 

[4o; 

[41 

[42; 
[43; 
[44; 
[45; 



Togan, E., Chu, Y., Imamoglu, A. & Lukin, M. Laser cooling and real-time measurement of the 

nuclear spin environment of a solid-state qubit. Nature 478, 497-501 (2011). 
Wilson- Rae, L, Zoller, P. & Imamoglu, A. Laser cooling of a nanomechanical resonator mode to 

its quantum ground state. Phys. Rev. Lett. 92, 75507 (2004). 
Mintert, F. & Wunderlich, C. Ion-trap quantum logic using long-wavelength radiation. Phys. 

Rev. Lett. 87, 257904 (2001). 
Cirac, J., Blatt, R., Parkins, A. & Zoller, P. Spectrum of resonance fluorescence from a single 

trapped ion. Phys. Rev. A 48, 2169 (1993). 
S0rensen, A. & M0lmer, K. Quantum computation with ions in thermal motion. Phys. Rev. Lett. 

82, 1971-1974 (1999). 
Bermudez, A., Schmidt, P. O., Plenio, M. B. & Retzker, A. Robust trapped- ion quantum logic 

gates by continuous dynamical decoupling. Phys. Rev. A 85, 040302 (2012). 
Lemmer, A., Bermudez, A. & Plenio, M. B. Driven geometric phase gates with trapped ions. 

arXiv:1303.5770 (2013). 
Roos, C. Ion trap quantum gates with amplitude- modulated laser beams. New J. Phys. 10, 

013002 (2008). 
Togan, E. et at. Quantum entanglement between an optical photon and a solid-state spin qubit. 

Nature 466, 730-734 (2010). 
Greffet, J.- J. et al. Diamond particles as nanoantennas for nitrogen- vacancy color centers. 

arXiv:1107.0502 (2011). 
Beveratos, A., Brouri, R., Gacoin, T., Poizat, J. -P. & Grangier, P. Nonclassical radiation from 

diamond nanocrystals. Phys. Rev. A 64, 061802 (2001). 
Neumann, P. et al. Excited-state spectroscopy of single NV defects in diamond using optically 

detected magnetic resonance. New J. Phys. 11, 013017 (2009). 
Lukin, M. & Hemmer, P. Quantum entanglement via optical control of atom-atom interactions. 

Phys. Rev. Lett. 84, 2818-2821 (2000). 
Bennett, S. et al. Phonon-induced spin- spin interactions in diamond nanostructures: application 

to spin squeezing. arXiv:130L2968 (2013). 
Hettich, C. et al. Nanometer resolution and coherent optical dipole coupling of two individual 

molecules. Science 298, 385-389 (2002). 
Waldermann, F. et al. Creating diamond color centers for quantum optical applications. Diam. 

Relat. Mater 16, 1887-1895 (2007). 
Tamarat, P. et al. Stark shift control of single optical centers in diamond. Phys. Rev. Lett. 97, 

83002 (2006). 
Santori, C. et al. Nanophotonics for quantum optics using nitrogen- vacancy centers in diamond. 

Nanotechnology 21, 274008 (2010). 
Neumann, P. et al. Quantum register based on coupled electron spins in a room-temperature solid. 

Nature Phys. 6, 249-253 (2010). 
Cohen- Tannoudji, C, Dupont-Roc, J. & Grynberg, G. Atom-Photon Interactions (Wiley- VCH, 

2004). 
James, D. & Jerke, J. Effective Hamiltonian theory and its applications in quantum information. 

Can. J. Phys. 85, 625-632 (2007). 



